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ABSTRACT 

Measurements of the total amount of stars locked up in galaxies as a function of host halo mass 
contain key clues about the efficiency of processes that regulate star formation. We derive the total 
stellar mass fraction fi, as a function of halo mass M^qqc from z = 0.2toz = l using two comple- 
mentary methods. First, we derive using a statistical Halo Occupation Distribution model jointly 
constrained by data from lensing, clustering, and the stellar mass function. This method enables us 
to probe f^, over a much wider halo mass range than with group or cluster catalogs. Second, we derive 
/* at group scales using a COSMOS X-ray group catalog and we show that the two methods agree to 
within 30%. We quantify the systematic uncertainty on f^, using abundance matching methods and we 
show that the statistical uncertainty on 10%) is dwarfed by systematic uncertainties associated 

with stellar mass measurements 45% excluding IMF uncertainties). Assuming a Chabrier IMF, we 
find 0.012 < < 0.025 at Afgooc = 10" Mq and 0.0057 < /* < 0.015 at A/sooc = 10^^^ M©. These 
values are significantly lower than previously published estimates. We investigate the cause of this 
difference and find that previous work has overestimated fi, due to a combination of inaccurate stellar 
mass estimators and/or because they have assumed that all galaxies in groups are early type galaxies 
with a constant mass-to-light ratio. Contrary to previous claims, our results suggest that the mean 
value of is always significantly lower than /gas for halos above 10^^ Mq. Combining our results with 
recently published gas mas fractions, we find a shortfall in /*-l-/gas at i?5ooc compared to the cosmic 
mean. This shortfall varies with halo mass and becomes larger towards lower halos masses. 
Subject headings: cosmological parameters, cosmology: observations, diffuse radiation, galaxies: clus- 
ters: general, galaxies: stellar content. X-rays: galaxies: clusters 



1. INTRODUCTION 

Dark matter halos contain baryons in the forms of stars 
and gas, but there is active debate about the level of 
agreement between the baryon fraction in halos and the 
cosmic mean, /f, = , as well as which baryonic phase 
dominates as a function of halo mass Mh- Cosmologi- 
cal simulations suggest that the baryon content of mas- 
sive galaxy clusters (M/i ^ 10^^ M©) should be within 
10% of the unive rsal value {e.g., IKravtsov et al.l [20051 : 
lEttori et al.l 120061 ). In galaxy clusters, most baryons 
reside in the hot, diffuse. X-ray emitting intra-cluster 
medium (ICM). Detailed X-ray measurements suggest 
that the hot gas fraction (/gas) is considerably lower 
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than the cosmic value in the inner regions of clusters 
(e.g.. lEttor i' 2003": Lin ct al.' 2003: Vikhhnin et al. '20061 
McCarthv et al 2007; A rnaud et al. 2007; Allen et aL| 
2nnilS un et al.H2009l ) though recent measurements sug- 
gest that /gas might approach th e universal value at 
larger radii (jSimionescu et aLll2011[ ). 

Depleted gas mass fractions at small radii might indi- 
cate that the ICM has been affected by non gravitational 
processes such as star formation or feedback from Active 
Galatic Nuclei (AGN). Indeed, sources of heat such as 
these that are located near the centers of clusters might 
pump thermal energy into the ICM causing it to ex- 
pand towards the outskirts (thus lowering /gas in cluster 
cores). Significant sources of heat might even be capable 
of rem oving gas from lower-mass halos (|McCarthv et al.l 
1201 It ). On the other hand, the fact that /gas is lower 
than expected at small radii might also simply indicate 
that some gas has been transformed into other baryonic 
components such as galaxy stellar mass or Intra-Cluster 
Light (ICL). Discriminating between these two scenarios 
is key towards understanding the efficiency of processes 
that regulate star formation and for tracing the thermo- 
dynamic history of ICL. 

Galaxy groups are highly interesting in this respect 
since they have shallower potential wells than clus- 
ters and should therefore be more sensitive to non - 
gravitational processes {e.g., IMcCarthv et"an [20T1 . 
Measurements of the baryon fraction at group scales, 
however, are disparate and subj ect to much debate (e.g., 
Lin et al.ll20 03: Gonzalez ct al. 20071 iGiodini et al.ll200a 
Balogh et al.l 120071 ). In particular, one unresolved issue 
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is whether or not the baryon content of galaxy groups is 
dominated by hot gas (hke clusters) or if other baryonic 
phases such as galaxy stellar mass or ICL also make a sig- 
nificant contribution. From X-ray studies there is a fairly 
clear consens us that /yas decrease s towards lower halo 
masses (e o.. iVikhlinin et al 1120061 fArnaud ct aL 120071 : 
iSun et al.H2009l: iPratt et al.ll2009D . The two other com- 
ponents, fi, and /icL, have proved more challenging to 
measure and are largely responsible for disagreements 
in the literature regarding the group scale baryon frac- 
tion. The ICL component is diffuse and faint (with a 
typical surface brightness of 27-32 mag arcsec"^ in the r- 
band) and therefore ex tremely difficult to measur e (e.g., 
Feldmeier et~all 120041 : [Lin & Mohr 2004; Zibctt i et alJ 
20051: IGonzalez et al.l [2b05i: iKrick et al.ii200a . The 
component is subject to uncertainties regarding the as- 
signment of galaxies to groups and clusters, but also 
to uncertainties inherent to stellar mass measurements 
themselves. 

Previous results on this topic can be broadly divided 
into two categories. In the first category, the decrease 
in /gas towards lower halo masses is found to be ex- 
actly compensated by an increased contribution from 
/*+/iCL in such a way that the baryon fraction is con- 
stant with hal o mass and is found to be either slightly 
below ( Gonza lez et al.l 120071 ) or significantly below ft 
(|Aiidreonii2010l ). This would imply that there is a simple 
trade-off between /gas and /*-I-/icl from group to clus- 
ter scales. In the second category, /*-f/icL do not ex- 
actly compensate for the decrease of /gas at group scales 
and thus the baryon fraction is found t o be mildly de- 
creasing towards lo wer halo masses fe.q- ILin et 11112001 
IGiodini et al.|[2009() . This could imply for example that 
the ICM is more strongly affected by feedback mecha- 
nisms in groups than in clusters. 

Most of the work to date has not directly probed fi,, 
but has focused instead on the relationship between halo 
mass and the total K-band luminosity of group and clus- 
ter galaxies (e.g., Ilin et al.l 120031: IRamella etaLl 120041: 
iBalogh et all 120071 7 1201011 Although K-band luminos- 
ity is most sensitive to the low-mass stars that dom- 
inate the total mass of most stellar populations, it is 
not a direct tracer of this mass. Other work, for lack 
of multiband photometry, has estimated fi, using simple 
mass-to- light ratio (M/L) estimates (e.g. iLin ct al. 2003' 
Gonzalez ct al. 2007; Balogh ct al. 2007; Lagan a et alJ 
2008 ; Giodiniet al. 2009; Andreon 2010; Balogh et al.l 
20ld iDai et al.ll2010l:"lLagana et al.H201lf ). In practice, 
however, mass-to-light ratio values can vary strongly 
with galaxy color, metallicity, and age. No single mass- 
to-light ratio can capture the complexity of the full 
galaxy population (Ilbert et al. 2009) . No studies to date 
have actually accounted for the full morphological mix of 
galaxies in group environments and used Spectral Energy 
Distribution (SED) fitting methods to derived /*. 

In this paper we measure /* using stellar masses de- 
rived from full SED fitting to multi-band photometry (in- 
cluding K-band) from the COSMOS survey. We tackle 
the issue of measuring fi, using two complementary ap- 
proaches. First, we use a Halo Occupation Distribution 
(HOD) model that has been calibrated to fit measure- 
ments of galaxy clustering, galaxy-galaxy lensing, and 
the galaxy SMF as a function of redshift. Using this 



method we can derive /* down to much lower halo masses 
(a few 10"^^ M0) than possible using group and cluster 
catalogs. Second, we use a group membership catalog to 
directly calculate /* for 10^^ < Mh/Mq < lO" by sum- 
ming together the stellar masses of group members. In- 
terestingly, we find results that are significantly different 
than previous work. Assuming a Chabrier Initial Mass 
Function (IMF) , our results suggest that previously pub- 
lished estimates of /,t on group scales are overestimated 
by a factor of about 2 to 5. The discrepancy is only 
partially reduced by assuming a Salpeter IMF. We in- 
vestigate the cause of this discrepancy and find that pre- 
vious work has overestimated /* due to a combination of 
inaccurate stellar mass estimators and/or because they 
have assumed that all galaxies in groups and clusters are 
passive early type galaxies with a constant mass-to-light 
ratio. 

The layout of this paper is as follows. The data are 
briefly described in Section[2]followed by the presentation 
of our method used to derive fi, in Section [S] Our main 
results are presented in Section|4l Finally, we discuss the 
results and draw up our conclusions in Section [S] 

We assume a WMAP5 ACDM cosmology with fin, = 
0.258, r^A = 0.742, n^h"^ ^ 0.02273 = 0.963 
0-8 = 0.796, Ho = 72 km s^^ Mpc'^ (jHinshaw et al.l 
[20091) . We assume a WMAP5 cosmo logy to maintain 
consistency with the HOD analysis of iLeauthaud et all 
(2011b), which is an key ingredient for this paper. All 
distances are expressed in physical Mpc units. Halo mass 
is denoted Mh, and we define M500C = M{< rsooc) = 
500pcrit|T5ooc where rsooc is the radius at which the 
mean interior density is equal to 500 times the critical 
density (pcrit)- Stellar mass, denoted M^,, is derived 
using a Chabrier Initial Mass Function (IMF) unless 
otherwise specified. The integrated (total) stellar con- 
tent is denoted /^ and is defined as /* = M'°*/M5ooc- 
Herein, /^ is the sum of two components; a contribu- 
tion from sateUite galaxies denoted fl^*- = Ml^^-fM^ooc 
and a contribution from central galaxies denoted as 
jcen = M;^™/M5ooc; it does not include the ICL. The 
function E{z) = H{z)/Hq = ^Jf^^^JT^^Yf'^^TlA repre- 
sents the Hubble parameter evolution for a flat metric. 
Stellar mass scales as 1/Hq. Halo mass scales as I/Hq. 
All magnitudes are given on the AB system. 

2. DATA DESCRIPTION 

The COSMOS survey fScoviUe et all[2007l ). centered at 
10:00:28.6, +02:12:21.0 (J2000), brings together a broad 
array of panchromatic observations with imaging data 
from X-ray to radio wavelengths and a large spectro- 
scopic follow- up program (zCOSMOS) with the VLT 
(jLillv et al.l 120071 ). In the following sections we describe 
the key COSMOS data sets and previously published 
analyses relevant to this paper. 

2.1. Stellar Mass Estimates 

The goal of this paper is to estimate the total stel- 
lar content of dark matter halos. Stellar mass esti- 
mates are therefore an important component of this 
paper. We adop t the same stellar masses as used by 
ILeauthaud et al.l (|2011bD (hereafter "Lll") and George 
et al. (2011) (hereafter, "Gel l"). These ha ve been de- 
rived in a similar manner as iBundv et all (201Q) but 
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are based on updated redshift information (vl.7 of the 
COSMOS photo-z catalog and the latest available spec- 
troscopic redshifts as compiled by the COSMOS team) 
and use a slightly different cosmology {Hq = 72 km 
s~i Mpc~i instead of = 70 km s ~^ Mpc^^). We 
refer to Lll and iBundv et al.l ()2010[ ) for details re- 
garding the derivation of the stellar masses and only 
give a brief description here. Stellar masses are calcu- 
lated using COSMOS ground-based photometry (filters 
u*, Bj, Vj, ,r~^ , z'^ , Kg) and the depth in all bands 
reaches at least 25th magnitude (AB) with the Ks-ha.nd 
limited to Kg < 24. Stellar ma sses are estimated u sing 
the Bayesian code described in iBundy et al.l (|2006[ ) as- 
suming a Chabrier In itial Mass Function (IMF) and the 
iCharlot fc Fal (|2000') dust model. The stellar mass com- 
pleteness limits are shown in Figure 2 of Lll and vary 
from lO^-S Mq at z 0.37 to lO^-^ Mq at z = 0.88. 

2.2. Cosmos X-ray group membership catalog 

The entire COSMOS region has been mapped through 
54 overlapping XMM-Newton pointings and additional 
Chandra observations cover the central region (0.9 
degrees^) to higher r esolution ( i Hasinger et al.l 120071 : 



iCappelluti eTaLlT2551 lElvis et al.n2009[ ). These X-ray 
data have been used to construct a COSMOS X-ray 
group catalog that contains 211 extended X-ray sources 
over 1.64 degrees^, spanning the range < z < 1, and 
with a rest-frame 0.1-2.4 keV luminosity range between 
IQ'*^ and 10'*'* erg s ~*. The general d a ta red uction pro- 
cess can be found in lFinoguenov et all ()2007D and details 
regarding improvement s and modifications to t he origi- 
nal catalog are given in lLeauthaud et ahl ()2010t ). George 
et al. (2011), and Finoguenov in prep. Hal o masses have 
been measured by iLeauthaud et al.l ()2010f ) for this sam- 
ple by using weak lensing to calibrate the relationship 
between X-ray luminosity (Lx) and halo mass. Groups 
in this catalog have halo masses that span the range 
10*3 < Mh/Mq < 10**. 

To estimate f^, from this group catalog, we need a 
method to distinguish galaxies that belong to groups 
from those in the field. For this we use the group mem- 
bership probability catalog from Gell. In Cell, the full 
Probability Distributio n Function (PDF) from the COS- 
MOS photo-z catalog (lllbert et al.l 120 091) has been used 
to derive a group membership probability for all galax- 
ies with F8UW < 24.2 from the COS MOS ACS galaxy 
catalog ()Leauthaud et al.l [20071 12011 bl The complete- 
ness and purity of the group membership assignment was 
studied in detail by Gell using spectroscopic redshifts 
and mocks catalogs. Within a radius of 0.5i?20Oc, the 
group membership catalog has an estimated purity of 
84% and a completeness of 92%. The completeness and 
purity of this catalog does not vary strongly with redshift 
at z < 1. 

2.3. Cosmos analysis of galaxy-galaxy lensing, 
clustering, and stellar mass functions 

In Lll a HOD model was used to probe the relationship 
between halo mass and galaxy stellar mass from z = 0.2 
to z = 1. Constraints were obtained by performing a 
joint fit (in three redshift bins) to galaxy-galaxy lens- 
ing, galaxy clustering, and the galaxy stellar mass func- 
tion, with a 10 parameter HOD mode L A detailed de- 
scription of this model can be found in ILeauthaud et al.l 



(|2011al) . In particular, as described in Section 2.3 of that 
paper, this model can be used to derive as a function 
of Mh- Hereafter, we refer to this approach as the "HOD 
method" . This method is discussed in more detail in 
Section O 

3. METHOD 

The aim of this paper is to derive fi, as a function of 
Mfi- To achieve this goal, we use two different and com- 
plementary measurements. First, we use a HOD model 
that has been calibrated to fit measurements of galaxy 
clustering, galaxy-galaxy lensing, and the galaxy SMF as 
a function of redshift. Second, we use a group member- 
ship catalog to directly calculate /* by summing together 
the stellar masses of group members. There are advan- 
tages and disadvantages inherent to each method. The 
HOD method might be considered to be a more indirect 
probe of however, it allows us to derive over a 
much wider halo mass range than with group and clus- 
ter catalogs. For example, in COSMOS, this technique 
allows us to probe down to the stellar mass complete- 
ness limit of the survey (M* ~ 7 x 10^ Mq at z = 0.37). 
The group catalog method is more direct and less prone 
to modeling uncertainties (for example, the assumption 
that satellite galaxies follow the same radial distributions 
as the dark matter). However, this technique is limited 
by the halo mass range for which group and cluster cat- 
alogs can be constructed {Mh > lO*'^ Mq) and is sub- 
ject to other uncertainties such as the accurate identi- 
fication of the centers of groups and cluster of galaxies. 
The group catalog method is also subject to uncertainties 
due to the completeness and purity of the membership 
assignment and to projection effects, that must be under- 
stood through mock catalogs. Finally, we compare both 
methods with the results fr om the abundanc e matching 
technique described in Beh roozi et al.l ()2010[ ). with the 
assumption of the COSMOS stellar mass function. This 
model has substantially fewer parameters but makes a 
specific assumption about the connection between galax- 
ies and dark matter subhalos. 

3.1. HOD method 

The relationship between halo mass and galaxy stel- 
lar mass can be described via a statistical model of 
the probability distribution P{N\Mh) that a halo of 
mass Mh is host to N galaxies above some thresh- 
old in luminosity or stellar-mass. This statistical 
model, commonly known as the halo occupation dis- 
tribution (HOD), has been considerably successful at 
interpreting the clustering properties of galaxies (e.g., 
iSeliak '2000; 'Peacock & Smith ^00^; 'Scoccimarro et al 
2001; Berhnd & Wcinbcrs 2002; Bullock et al. 2002 
Zeha vi et al.l j2002, 2005; Zhcne ct al. 120051 12007 



i Tinker et al.l[2007: Wake et al., 2011: Zehavi et al.ll2011 
IWhite et alTboilD . Although the HOD is usually in- 
ferred observationally by modeling measurements of the 
two-point correlation function of galaxies, £,gg(r), it can 
also be used to model other observables such as galaxy- 
galaxy lensing and the stellar- mass function. In Lll the 
HOD model was used to probe the relationship between 
halo mass and galaxy stellar mass from z = 0.2 to z = 1. 
Constraints were obtained by performing a joint fit to 
galaxy-galaxy lensing, galaxy clustering, and the galaxy 
stellar mass function. 
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In this paper, we adopt the parameter fits from Table 
5 in Lll and we use these to study as a function of 
Mh. We calcul ate fj, following t h e proc edure outlined in 
Section 2.3 of iLeauthaud et all (|2011al ). One potential 
concern with this method is that the calculation of /^^* 
might require extrapolating the HOD model beyond the 
lower and upper stellar mass bounds for which the model 
has been calibrated. As shown in Lll, f^, is dominated 
by the stellar mass of central galaxies for Mh < 10^^ 
Mq . This concern is therefore only valid for halos larger 
than about 10^'^ M0 where satellite galaxies represent 
the dominant contribution to fi,. 

We now investigate the stellar mass range of satellite 
galaxies that build up the bulk of fl^'^ (the satellite con- 
tribution to /*) in halos above Mh > lO" Mq. The 
original HOD analysis of Lll used M2006 halo masses 
(defined with respect to 200 times the mean matter den- 
sity p). For simplicity, for the tests in this paragraph 
we will also use M200& but the conclusions are also valid 
for M500C masses. Let us consider the stellar fraction for 
satellites in a fixed stellar mass bin, /f''*(A//i|M", M"), 
where M*^ represents a lower stellar mass limit and Af*^ 
represents an upper stellar mass limit. The expression 
for fT\Mh\Mt\ Mf") is given by: 

1 r*^" 

/r*(M,,|Mf ,Mf) = — / ^siM,\Mh)M,dM,. 

(1) 

where represents the condition al stellar mass 

functio n for satellite galaxies {e.g., ILeauthaud et al.l 
I2011al) . Using this expression, we have tested how 
fl'^^{Mh\Ml^,Mf) varies with the integral limits 
and M*^. We find that at fixed halo mass, the bulk of 
fl^'^ comes from galaxies with > 10^ M©. For exam- 
ple, for Maoob > 10^^ M©, 95% and 90% of /f * is due to 
satellites with Af » > 2 x 10^ M© and Af, > 7 x 10^ M© 
respectively. The HOD model of Lll is calibrated down 
to 6 X 10*^ M© at z = 0.37 and to 6 x 10^ M© at 2 = 0.88. 
Therefore, the bulk of fl^'^ arises from satellites that are 
well within the test ed limits of our mode l (for a similar 
calculation, also see iPuchwein et~aLll2010t) . 

As mentioned previously, the Lll analysis assumes 
A^aoofc halo masses whereas in this paper we use Afsooc 
halo masses so as to facilitate comparisons with previ- 
ous work. We therefore need to convert the Lll re- 
sults from Af2006 to Afsooc- For this we make two 
assumptions. First, we assume that satellites galax- 
ies are distributed according to the sa me profiles as 
their parent halos, namely, NFW profiles (Navarro et al.l 

19971) with g. t = Ch.io (iNagai fc K ravtsov 2005; 

Wetzel fc Whit^ MWl ITinker et aLlTfouD . Second, we 
assume that the distribution of satellite galaxies is inde- 
pendent of stellar mass (i.e, there is no "satellite mass 
segregation"). This assumption is supported by multi- 
ple observations jP racv et al. 2005; Hudson et al.l [2OIOI : 
Ivon de r Linden et al. 2010; We tzel et al.i i201ll ) but see 
also [van d en Bosch et al . (2008t)- Under these two as- 
sumptions, Z^'^* is unchanged by any halo mass con- 
version (because the same conversion factor applies to 
both the numerator and the denominator in fl'^^). The 
Z^®'^ component does vary, however, because the halo 
mass (denominator) changes while the mass of the cen- 



tral galaxy (numerator) remains constant. For details 
on how to convert hal o masses assuming a N FW pro- 
file, see Appendix C in lHu fc Kravtsov! (|2003D . For this 
conversion, we assume the m ass-concentration relation of 
iMufioz-Cuartas et al.l ()2011h . 

3.2. Group catalog method 

In this section, we describe our procedure for calculat- 
ing fi, at group scales using the COSMOS X-ray group 
catalog. Figure [1] shows the group sample as a function 
of redshift and Lx • Blue boxes show the binning scheme 
that we adopt for this paper. The redshift limits are se- 
lected to match those of Lll and are zi — [0.22,0.48], 
Z2 = [0.48,0.74], and 23 = [0.74, 1]. 

To ensure a robust group sample with a relatively clean 
member list, we include only groups with reliable opti- 
cal associations having greater than 3 members, and we 
exclude close neighboring systems and groups near the 
edges of the field or masked regions where member as- 
signment is difficult (catalog flags XFLAG=1 or 2, and 
POOR=MERGER=MASK=0; See Cell for flag definitions). 
This ensemble of quality cuts is synthesized by a global 
flag in Cell called flag .include (we select groups with 
FLAG_INCLUDE=1). In total, after these cuts, our catalog 
contains 129 groups at z < 1. 
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Fig. 1. — Cosmos X-ray group sample as a function of red- 
shift and Lx E{z)~^ , where Lx represents X-ray luminosity and 
E(z) = H{z)/Ho. Large blue circles indicate groups included in 
this analysis, and small black circles show groups that are excluded 
from this analysis (see text for selection flags). Yellow diamonds 
show systems where the centering is deemed uncertain based on 
the analysis of George et al. in prep., because the most massive 
group member lies far from the X-ray center. Blue boxes indicate 
the binning scheme. The redshift limits for this binning scheme 
are chosen to match the HOD analysis of Lll. 

The determination of group centers is an important 
and non-trivial task, especially since these structures do 
not always have a visually obvious central galaxy. Mis- 
centering effects could lead to an underestimate of f^,, es- 
pecially if the mis-centering is large enough that the true 
central galaxy is excluded from the group selection for 
example. George et al. in prep show that weak lensing 
is an excellent tool that can be used to optimize center- 
ing algorithms by maximizing the stacked weak lensing 
signal at small radii (below about 500 kpc). For exam- 
ple, George et al. in prep also show that various naive 
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centering algorithms such as the centroid of galaxy mem- 
bers, even if weighted by stellar mass or luminosity, is a 
poor tracer of the centers of dark matter halos. We use 
the group centers from George et al. in prep that are 
found to maximize the stacked weak lensing signals at 
small radii, namely, the most massive group galaxy lo- 
cated within an NFW scale radius of the X-ray centroid, 
denoted MMGGscaie- We note, however, that George et 
al. in prep have also identified a sub-set of groups (of 
order 20%) that host a more massive galaxy at a greater 
distance from the X-ray centroid than the NFW scale ra- 
dius. For these groups, the weak lensing does not clearly 
indicate a preferred center among the two options. We 
include these systems in our analysis but we have tested 
that our results are not affected by this choice. 

Group members are selected according to the following 
criteria: 

1- ^incm > 0.5 where Pmom denotes the group mem- 
bership probability derived in Gell. 

2. r2d < ^500c where r2d is the projected distance to 
the group center (defined hereafter by the position 
of the MMGGscaic). 

3. Galaxies must have a stellar mass above the com- 
pleteness limit of the bin in consideration (see Table 
|T]). Stellar mass completeness limits for the group 
membership catalog are higher than those quoted 
in Section [2TT1 because the group membership cat- 
alog is limited to F8UW < 24.2. 

Group members are selected only within a radius of 
-R500C for two reasons. Gell have shown that the pu- 
rity of the group membership assignment drops at larger 
radii (see their Figure 5). By only considering galaxies 
at r2d < -R500C, we reduce the errors due to purity cor- 
rections. In addition, most previously published work on 
this topic has used r2d < Rbooc] this choice thus facili- 
tates comparisons with previous work (see i 34.5|) . 

Once group members have been selected following the 
criteria outlined above, fi, is derived for each group as the 
sum of the stellar masses of all group members associated 
with this group divided by the halo mass. Halo masses 
are derived from the X-ra y luminosity following th e weak 
lensing scaling relation of iLeauthaud et ahl (|201Clf ). 

Although in principle more direct than the HOD 
method, measuring estimates from the group catalog 
has systematic errors that must be corrected for. The 
two main effects on the measurement are: 

1. The completeness and purity of the membership 
selection. Here we include both contamination in 
the group membership selection due to neighbor- 
ing galaxies but also the de-projection of a spher- 
ical NFW profile. Indeed, the quantity that we 
are interested in is contained within a sphere of 
radius -Rsooc whereas our membership selection is 
akin to selecting galaxies in a cylinder in redshift 
space. We therefore include a correction factor for 
this de-projection. 

2. The stellar mass completeness of the group mem- 
bership catalog (i.e., we need to account for the 
contribution to f^, from galaxies below the com- 
pleteness limit). 



To derive these correction factors, we use a suite of 
COSMOS-like mock catalogs that have been described 
in Gell. Mocks were created from a single simulation 
(named "Consuelo") 420 Mpc on a sideQ, resolved 
with 1400^^ particles, and a particle mass of 1.87x10® 
h^^ Mq. This simulation can robustly resolve halos with 
masses above ^ 10^^ Mq which corresponds to a 
central galaxy stellar masses of ~ 10*'^ Mq, well- 
matched to our completeness limit of F814W^ = 24.2 at 
low redshift. This simulation is part of the Las Damas 
suitcQ (McBride et al. in prep). Halos are identified 
within the simulation using a friends-of-friends finder 
with a linking length oi b — 0.2, and populated with 
mock galaxies according to the HOD model of Lll. See 
Gell for further details regarding these mock catalogs. 

We derive a single overall correction factor to due to 
the two effects described above. To do so, we apply the 
group membership selection to the mock catalog and we 
compare the measured value of (in the mocks) to the 
true value of fi, (within a sphere of radius i?5ooc)- We 
derive correction factors separately for the contribution 
to from satellite and central galaxies. These correction 
factors are denoted Csat and Cccn respectively. Csat is 
defined as: 

Csat = logio(Mf °^'^-*™'VAC°'''"'=') (2) 

where j\,/™ock.truth represents the "true" total mass from 
the mocks in a sphere of i?5ooc and down to the com- 
pleteness limit of the mocks and jv^mockjsci represents our 
group member selection applied to the mocks. Ccon rep- 
resents the equivalent for central galaxies. The correction 
factors are listed in Table[2]for the three redshift bins and 
are at most 14%. The completeness and purity correc- 
tion component of Csat is negative (the measured value 
of /^^' is higher than it should be because of contam- 
ination) whereas the stellar mass completeness compo- 
nent of Csat is positive (the measured value of f^'^^ must 
be augmented to account for satellites below our com- 
pleteness limit). The magnitude of the two corrections 
is fairly similar, but the opposite signs mean that these 
corrections tend to cancel each other out. The redshift 
trends cause the net correction for Csat to be negative at 
z ~ 0.37 and positive at z ^ 0.88. 

Due to the finite mass resolution of our mock catalogs, 
the contribution to f^, from very low mass galaxies (those 
that arc not contained in our mock catalogs) will not 
be included. Using Equation [TJ we estimate that this 
contribution should only be of order 1-2%, and we thus 
neglect it here. 

4. RESULTS 

4.1. Results from the HOD method 

We now calculate f^, using the best fit HOD parame- 
ters from Lll for each of the three redshift bins. The 
results are shown in Figure [H The solid dark blue, blue, 
and turquoise lines show as a function of M500C for 

^ In this paragraph, numbers are quoted for Hq = 100 h 
km s~^ Mpc~^. The assumed cosmology for Consuelo is Qm = 
0.25, r^A = 0.75, n^h^ = 0.02273, = 1.0, o-g = 0.8, Hq = 70 
km Mpc~l. 

^ D etails regarding this simulation can be found at 
|http : //Iss . phy . vanderbilt . edu/lasdcunas/simulat ions .html | 
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TABLE 1 

Group binning scheme 



Bin ID Redshift lower hx limit upper Lx limit M« completeness limif 

logio(LxE(z)-Verg s-i) logio(LxE(z)-Verg s'l) logio(M./M0) 



A 


21 = 


0.22,0.48 


41.96 


42.43 


9.3 


B 


Zl = 


0.22,0.48 


42.43 


43.2 


9.3 


C 


Z2 = 


0.48,0.74 


42.43 


43.2 


9.9 


D 


23 = 


0.74, 1] 


42.76 


43.5 


10.3 



The stellar mass completeness limits cited here are for the group membership catalog. 



TABLE 2 

Correction factors for group catalog 



Bin ID 


Redshift 


C'sat 


C'ccii 


A 


21 = [0.22,0.48] 


-0.029 


-0.019 


B 


21 = [0.22,0.48] 


-0.066 


-0.004 


C 


22 = [0.48,0.74] 


-0.047 


-0.017 


D 


23 = [0.74, 1] 


0.003 


-0.017 



z ~ 0.37, z ^ 0.66, and z ~ 0.88 respectively. The con- 
tribution to fi, from central galaxies, f^°'^, is shown by 
the dotted red line and the contribution from satellites, 
Z^*^*, is shown by the dashed yellow line, fi, is dominated 
by central galaxies at Msooc ^ lO^'^'^ Mq and by satel- 
lites at M500C > 10^^-^ Mq. Note that /* in this figure 
does not include ICL. 

The redshift evolution of has been discussed at 
length in Lll. Briefly, the "pivot halo mass" is de- 
fined as the halo mass where reaches a maximum 
(Afsooc — 10^^ M0). At fixed halo mass and below the 
pivot halo mass scale, increases at later epochs. This 
evolution can be explained by the fact that in this regime, 
central galaxy growth outpaces halo growth. At fixed 
halo mass and above the pivot halo mass scale, there 
is tentative evidence that /* declines at later epochs. 
In Lll, we suggest that this trend might be linked to 
the smooth accretion of dark matter, which brings no 
new stellar mass, and amounts to as much as 40% of the 
growth of dark matter halos fe.g.. iFakhouri &: ]V[all2010[ ) 
and/or the destruction of satellites and a growing ICL 
component. 

The shaded grey area in Figure [5] indicates the statis- 
tical error on /^,. The statistical error is small compared 
to the systematic uncertainty associated with the actual 
derivation of stellar masses. We will return to this issue 
shortly in Section 14.31 

4.2. Results from the group catalog 

We now calculate f^, at group scales using the COS- 
IVIOS X-ray group catalog. Figure |3] shows our estimates 
for /^°" (top row), J^'^' (middle row), and (bottom 
row) as derived from the group catalog (the correction 
factors from Table [2] have been applied) . Each colored 
data point in Figure [3] represents a measurement from 
one X-ray group. The thick black data point repre- 
sents the mean value in each halo mass bin (see Table 
[T]). The horizontal error corresponds to the width of the 
bin while the vertical error represents a / ^^Ng where a is 
the measured dispersion in the bin and Ng is the number 
of groups in the bin. 

The solid filled symbols in Figure |3] indicate systems 



for which the weak lensing analysis of George et al. in 
prep suggests that the group centering might be unreli- 
able (MMGGscaio is different than MMGGr2oo)- We do 
not find a strong correlation between these systems and 
either halo mass or /^,, and find that the results are un- 
changed if these systems are removed from the analysis. 

We find excellent agreement between the HOD method 
and the group catalog method for Z^''" (compare the red 
squares to the red dotted line). Regarding /^^', there 
is a small but consistent offset between the two meth- 
ods (compare the yellow diamonds to the yellow dashed 
line). Indeed, in all three redshift bins, the f^^^ values as 
measured from the group catalog are higher than those 
derived from the HOD method by about 20% to 50%. 
As mentioned in the previous section, the HOD predic- 
tions for /^^* in Figure [3] have been derived under the 
assumption that the distribution of satellite galaxies fol- 
lows the same NFW as their parent halos. We hypothe- 
size that the offset for /^^* seen in Figure |3] might indicate 
that either a) the satellite mass distribution has differ- 
ent concentration than the dark matter or b) the satellite 
mass distribution varies with stellar mass ("mass segre- 
gation"). The observed level of offset, however, does not 
affect the main conclusions of this paper and so we leave a 
more detailed investigation of this effect for future work. 

Overall, the two methods show a remarkable level of 
agreement especially considering the fact that they are 
subject to very different types of systematic errors. The 
two methods agree to within 30% for the estimated value 
for at group scales. 

4.3. Systematic errors associated with f^, measurements 

The statistical uncertainty on fi, (from both the HOD 
method or the group catalog) is much smaller than 
the systematic uncertainty associated with the deter- 
mination of ste l lar m asses. For example, according to 
iBehroozi et al.l ()2010D . the typical systematic error as- 
sociated with photometric-based stellar mass estimates 
are of order 0.25 dex. This systematic uncertainty arises 
from the choice of a stellar population synthesis (SPS) 
model, the choice of a dust attenuation model, and the 
assumed functional form of the star formation history. 
In addition, there are also further uncertainties due to 
the choice of an IMF. Further details regarding system- 
atic errors in stellar m ass measurements can be found in 
IBehroozi et all (poTol ) and IConrov eFall (pOOl . 

In this section, we attempt to quantify how the sys- 
tematic errors described above might translate into an 
uncertainty on f^,. Note that the aim here is only to 
provide a first attempt to discuss systematic errors on 
due to uncertainties in stellar mass estimates (to our 
knowledge, this has not yet been quantified). We will 
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Fig. 2. — Total stellar content in galaxies (not including ICL) as a function of halo mass. Results from the HOD method constrained 
by COSMOS data are compared to the cosmic baryon fraction measured by the Wilkinson Microwave Probe (WMAP5; = = 

0.171 ± 0.009; IDunklev et~all |2009|) . The results for three redshift ranges zi, 22, and 23 are shown by the solid dark blue, blue, and 
turquoise lines respectively. The shaded region shows statistical errors only. Systematic errors are quantified in Section 14.31 (also see Figure 
O. /* is dominated by central galaxies at Msooc ^ 10^^'^ Mq (dotted lines) and by satellites at Msooc lO"*^"^'^ Mq (dashed lines). 



make certain assumptions that will undoubtedly need to 
be refined in future work. We focus on the low redshift 
results here (2 < 0.4) because the method used here re- 
lies on the availability of a large number of published 
SMFs. 

To begin with, we note that the SMF provides the 
strongest observational constraints on the HOD model of 
Lll. The other two observables used by Lll (clustering 
and galaxy-galaxy lensing) contain additional informa- 
tion, but this is less dominant than the information from 
the SMF. In this paper, we will therefore make the sim- 
plifying assumption that the main source of systematic 
error in the determination of is due to uncertainties 
related to the SMF. Under this assumption, the task of 
determining systematic errors on fi, now becomes one of 
determining systematic errors on the observed SMF. We 
consider two distinct sources of error for the SMF: 

1. A systematic error due to the choice of an IMF. 
In this paper, we specifically consider the Salpeter 
and the Chabrier IMF (Salpeter 1955; ChabriiB 
I2OO31V To first order, the choice of a Chabrier ver- 
sus a Salpeter IMF, yields masses that are lower 
by about 0.25 de>|3. For simplicity, we will assume 
here that we can simply convert between these two 
IMFs using a difference of 0.25 dex and we will re- 
fer to this difference as the "IMF systematic error 
on the SMF" . Note that anoth er IMF that is often 
considered is the Kroupa IMF (|Kroupall200ll) . The 

^ In practice, this conversion depends on the adopted SPS model 
and also the star formation history. The conversion can vary at the 
0.05 dex level: logio(M^""V^F'""'') = 0.25 ± 0.05 dex. 



Kroupa IMF yields stellar masses that are larger 
than a Chabrier IMF by about 0.05 dex. We ne- 
glect any possible variations of the IMF with either 
galaxy type and/or redshift but note that this is an- 
other potential (and currently poorly determined) 
source of systematic error on f^. 

2. All other sources of systematic error on the SMF, 
for example, those due to the choice of a SPS 
model, the choice of a dust attenuation model, and 
the assumed functional form of the star formation 
history. We will refer to these together as the "non- 
IMF systematic errors on the SMF" . 

To estimate the non-IMF systematic errors, we con- 
sider a compilation of various published SMFs at 
z < 0.4. Figure |4] shows an ensem ble of low-z 
(z < O A) mass functio n s froni COSMOS (jProrv et al.l 
I2OO9I: iLeauthaud et al.' '201 IbD and from the SPSS 
(jPanteTet a l. 2007; Bald rTet al.l [20081 iLi fc White! 
|2009| ). All mass functions in Figure |4] have been con- 
verted to our assumed value of h = 0.72 and to a 
Chabrier IMF. Observed variations between these differ- 
ent mass functions are due to a combination of system- 
atic error, measurement erroi0, sample variance, and red- 
shift evolutior0. We will adopt a conservative approach 

A higher level of measurement error will lead to an inflated 
stellar mass function at the high mass end due to Eddington bias; 
see for example discussion in Bchroozi et al. 2010. 

^ Most studies flnd very little redshift evolution in the total SMF 
out to z = 1 (e.o.. lBundv et al.ll2006D . Therefore, although there is 
a considerable time span between 2 = and z = 0.4, the differences 
between COSMOS and SDSS SMFs over this time period are likely 
primarily due to systematic errors and not redshift evolution. 
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Fig. 3. — Comparison between inferred from our HOD model with results from the COSMOS X-ray group catalog as a function of 
redshift. The HOD values for and /* are shown respectively by the dotted red line, and the dashed yellow line, and solid 

blue line. The red squares in the three top panels show Z^*^" measurements from the COSMOS X-ray group catalog (each data point 
represents one X-ray group) for the three redshift bins. Black data points with errors show the mean value in each bin. The horizontal 



error corresponds to the width of the bin while the vertical error represents a/^ 



Here a is the measured dispersion in the data and Ng 



is the number of groups in the bin. Yellow diamonds in the three middle panels show /|^* , and blue diamonds in the lower three panels 
show f*=f^™+f^^^- Grey solid points indicate systems for which the centering is ambiguous (excluding these systems does not change our 
results) . 

and assume that all the observed variation in Figure |4] is 
due to non-IMF systematic errors in the determination 
of stellar masses (this is an assumption that will need 
to be improved on in future work). We define a lower 
envelope and an upper envelope that is designed to en- 
compass the observed range of SMFs. This is shown by 
the grey shaded region in Figure |4l In what follows, we 
assume that this grey shaded region represents the non- 
IMF systematic errors on the SMF. There is, however, 
one caveat with this method. Our non-IMF systematic 
error margin might be underestimated due to the fact 



that not all SPS models are represented by the set of 
SMFs shown in Figure |^. For example, the SPS mod- 
els of iMarastonl (|2005f ) are not represented in Figure 
To date, however, ther e are no publ ished z < 0.4 SMFs 
that have used the Ma rastonl ()2005[ ) models. We there- 
fore simply note that this is a limitation with our current 



IPanter et al.i 
els. IBaldrv et al 



estimates from 



, , . u se lBruzual fc CharlotI (2003) SPS mod- 
_ mm use the average of four stellar mass 
Kauffmann et al. (2003), Pantcr et al. {200'%^ 
(2004), and Gallazzi et al. (2005). Li & WliiS 
" ( 200Ij). iDrorv et ail 

1 20091 ) and this paper use IBruzuaOrCh arlot ( 20CT1 SPS models. 



IGlazebrook et al 

II2009I) use masses from iBlanton fc Roweis 
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Fig. 4. — Compilation of stellar mass functions at z < 0.4. Sys- 
tematic error on the SMF not including the IMF is indicated by 
the shaded grey region. Variation due to measurement error, sam- 
ple variance between surveys, and redshift evolution is implicitly 
included, but we do not explore the full range of SPS models and 
SPS model parameters. 



The next step is to understand how these non-IMF 
errors on the SMF might translate into a systematic 
error on At present, the HOD model of Lll can only 
be used in conjunction with clustering and galaxy-galaxy 
lensing measurement. For this reason, we do not use the 
Lll HOD method for this step. Instead, we will adopt 
an "abundance matching" methodology that assumes 
there is a monotonic correspondence between halo 
mass (or circular velocity) and galaxy stellar mass (or 
lumin o sity) (e.g.. iKravtsov e t al. 2004; Vale & Ostriker 



2004 iTasitsiomi et al.l 120041: I Vale k Ostrikeil 

ConroY et alj 120061: iConrov fc Wechsleii 

Moster et all l2010t iBehroozi et al.l 120101: iGuo et al.l 



20101 ). Abundance matching techniques are generally 



designed to fit the SMF alone and have fewer free 
parameters than the HOD meth od. We follow th e abun- 
dance matching techniques of IBehroozi et al.l ()2010[ ). 
using M500C halo masses, using both the upper and the 
lower limits of the grey shaded region in Figure |4l The 
abundance matching results are then reported in Figure 
[5] as the non- IMF systematic error margin for f^, (solid 
light grey region). The systematic errors margin for 
fi, assuming a Salpeter IMF are estimated by shifting 
this region upwards by 0.25 dcx in stellar mass (hashed 
region in Figure [5]). 

In Figure El the COSMOS resuhs are slightly higher 
than our non-IMF systematic error margin in Figure [5] at 
M5Q0C > 10^^ Mq. This is due to the fact that the HOD 
method of Lll and t he abundance matching method of 
IBehroozi et al.l (|2010D yield slightly different predictions 
for as function of halo mass; this is discussed further 
in the following section. 

4.4. Comparison w ith the abundance ma tching results of 
Wehroozi et~al\ (MM) 

IBehroozi et al.l ()2010D have derived constraints on 
(see their Figure 10) by abundance matching to the SDSS 
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Fig. 5. — Systematic error on due to uncertainties in stellar 
mass estimates. The dark blue line shows the COSMOS HOD re- 
sult at 2; ~ 0.37. The dotted blue line shows derived by applying 
the same abundance matching technique as iBehroozi et al.l l|2010l ) 
to the COSMOS z ~ 0.37 SMF. The HOD method and the abun- 
dance matching method yield slightly difference results at the 15% 
to 20% level. The dashed magenta line shows as measured by 

throozi et al. (2010) by abundance matching to the lLi fc Whitg 
)09) SDSS SMF. The difference between the CO SMOS HOD 
result (dark blue solid line) and the JBchroozi et al.| j201CI|) result 
(dashed magenta line) is largely driven by differences between the 
COSMOS mass function and the lLi fc~ White (2009) mass function 
(see Figure |4]l. The solid light grey region shows our the system- 
atic error margin for /*, assuming a Chabrier IMF. The fact that 
the COSMOS results lie above this region at A/500C > 10^'^ Mq is 
due to a 15% model uncertainty between the HOD method and the 
abundance matching method. The hashed region shows the sys- 
tematic error margin for /* assuming a Salpeter IMF (note that 
the hashed region simply corresponds to a 0.25 dex upward shift 
of the solid light grey region). 



SMF of lLi fc Whiti (|2009f) . Figure [5] shows these results, 
here translated to M500C halo masses (n iagenta dashed 
line) and assuming a Chabrier IMF. The IBehroozi et al.l 
([2010) results yield a lower amplitude for fi, than the Lll 
COSMOS results. This is largely driven b y differences 
betwe en the COSMOS mass function and the lLi fc White! 
(|200a mass function (shown in Figure |3]). 

In order to investigate how much of these differences 
are in fact due to differences in the assumed stellar mass 
function, we have applied the abundance matching tech- 
nique of Behroozi et al. (2010) to the COSMOS z - 0.37 
SMF. This is shown as the blue dotted line in Figure 
El This can then be directly compared with the Lll 
HOD results. We find that the HOD results in roughly 
15% higher values for all halo masses than the abundance 
matching technique shown here. 

There are several possible sources for the 15% discrep- 
ancy between the HOD and abundance matching meth- 
ods. One possibility is that the weak lensing and/or clus- 
tering is providing more information in the HOD analyses 
than is available given the stellar mass function alone, 
and that this information pushes the values slightly 
higher. Another potential issue is that conversions have 
been made in each case to use M500C halo masses (HOD 
results were calculated initially using M2oob and abun- 
dance matching was done using Myir); these conversions 
necessarily make assumptions about the radial profiles 
of satellites which may be inaccurate at the several per- 
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cent level. Another possibility is that the abundance 
matching method might suffer from satellite incomplete- 
ness in the N-body simulations (see e.g. Wu et al. in 
prep), however, this is unlikely to be important below 
the group mass scale. We leave tracking down the ex- 
act source of the difference to future work, and simply 
note that there is roughly a 15% level uncertainty in the 
determination of the non-IMF systematic uncertainty on 
fi, due to the method used to derive ft, from the SMF. 
Overa ll, the difference between Lll and iBehroozi et aD 
(j2010( ) is consistent with our systematic error margin for 
and these differences do not impact any of our con- 
clusions significantly. 

4.5. Comparison with previous results from group and 
cluster catalogs 

We now compare our f^, m easurements wi t h two 
recently published resu lts from iGonzalez et all (|2007[ ) 
(hereafter "Gz07") and iGiodini et al.l (|2009ir rhereafter 
"Gi09") which have been derived using group and clus- 
ter catalogs. Figure [6] shows our predictions for f^^ as 
a function of M^qoc compared to the measurements of 
Gz07 and Gi09. In both cases, our results suggest sig- 
nificantly lower values for f^, than reported by either of 
these previous studies. We now provide a more detailed 
comparison with both Gz07 and Gi09 and discuss possi- 
ble explanations for the source of the discrepancy. 



4.5.1. Comparison wit MGonzalez et all l200'i) 

The left panel of Figure [6] shows the comparison of 
our low-z results with Gz07. Yellow diamonds show esti- 
mates from Gz07 for the contribution to fi, from satellite 
galaxiesQ whereas dark blue squares show the total stel- 
lar fraction {fl'^^+ff'^'^+fY'^^)- The yellow diamonds in 
Figure |6] are therefore directly comparable to our predic- 
tion for fl^'^ which is shown by the yellow dash-dash line. 
Gz07 find that fl^'^ increases towards lower halo masses 
whereas we predict the opposite trend. It is clear from 
this figure that a large fraction of the observed discrep- 
ancy is due to the ff^ component (and not so much to 
jBCG jiCL-j^ ^-^Yl now investigate the source of 
this difference in further detail. 

To estimate stellar masses, Gz07 have used a Cousins 
I-band mass-to-light ratio {M/Lj = 3.6) calibrated from 
the dynamical m odeling of 2- d kinematic dat a from 
iCappellari enil ([1)06). The iCappellari erall (pOOel ) 
sample is comprised of early type galaxies from the 
SAURON sample (Bacon et al. 2001). We note two 
important caveats with this approach: 

1. The SAURON galaxy sample is comprised solely 
of elliptical (E) and lenticular (SO) galaxies which 
have larger M/L values on average than intermedi- 
ate and late type galaxies. If the early type fraction 
increases with halo mass as sugges t ed by several 
studies (e.g., [ Wcinma nn et al.l 120061 : iHansen et al.l 
120091 : iWetzel et al.. .20111 ). then this could explain 
both why Gz07 find a steeper slope and a higher 
amplitude than we do for f^^^ at group scales. 



Ttsat jg calculated from Table 1 (column 7) in IGonzalez et aLl 
II2007I ) using a Vega solar luminosity of Msun = 3.94 (A.Gonzalez, 



2. The M/ L ratios presented in ICappellari et al.l 
()2006D are dynamical (i.e total) and are there- 
fore sensitive to the dark matter fraction within 
the effective r adius (Re)- The dynam ical mass-to- 
light ratios in ICappellari et aD ()2006( ) might there- 
fore be biased high by up to 30% compared to 
the stellar mass-to-lig ht ratios (see Figure 17 in 
ICappellari et al.ir2006l for example). 

We now investigate the mass-to-light ratios of our 
group members compared to the value used by Gz07. 
For this exercise, we consider group members from the 
membership catalog of Cell in our two low redshift bins 
(bins A and B in Figure [T]). We derive Subaru z+-band 
mass-to-light ratios from our stellar mass catalog using 
the absolute i'^-band AB magnitude that is provided by 
the COSMOS photoz catalog. For this exercise we adopt 
an AB solar luminosity of Msnn = 4.5^3. Figure [7] shows 
the M/Li+ values for COSMOS group members. We use 
unextincted rest frame magnitudes from the COSMOS 
photoz catalog to divide the histogram in Figure [7] by 
galaxy color. More specifically, we define galaxy color as 
C = M{NUV) - m{r) where M{NUV) and M{R) are 
the unextincted rest-frame template magnit udes in the 
near-u ltraviolet and the R bands defined in jlbert et al.l 
(|2OT0l ). We adopt the same division as in that paper, 
namely: 

C < 1.2 "high activity" or "blue" 
1.2 < C < 3.5 "intermediate activity" or "intermediate" 
C > 3.5 "quiescent" or "red" 

We compute the difference in magnitudes between our 
Subaru i'^ filter and a Cousins I filter for a range of 
stellar population templates. We find that a Cousins 
M/Li = 3.6 corresponds to a Subaru i"*" mass-to-hght 
ratio of 3.67 < M/Li+ < 4.1 (the exact value depends 
on galaxy color). These values are represented in Figure 
[7| by the hashed magenta region. 

Figure [7] shows that the stellar masses inferred by us- 
ing a similar M /L as Gz07 are larger than our SED stel- 
lar masses by a factor of 2 to 10 (assuming a Chabrier 
IMF). The difference is particularly large for active and 
intermediate type galaxies. The large spread observed in 
Figure [7] highlights the fact that the group population is 
not well represented by any single ^"'"-band M / Li^ ratio 
value. Part of the observed spread can also be explained 
by the sensitivity of the i"'"-band luminosity to recent star 
formation, underscoring the need for NIR data to more 
accurately determine galaxy stellar mass. The differences 
between our mass-to-light ratios and the value used by 
Gz07 is reduced by assuming a Salpeter IMF (Figure [71 
right panel). Even with a Salpeter IMF, however, the 
Gz07 M/L ratio is still well above the mean M/L ratio 
of our galax y sample. This i s like ly to be due to the 
fact that the lCappellari et all ()2006l ) mass-to-light ratios 
are dynamical and are thus sensitive to dark matter in 
addition to stellar mass. 

Another important fact worth highlighting in Figure |7| 
is that a non negligible fraction of group members are 
star forming/intermediate type galaxies (also see Gell 



priv. comm) and a Cousins I-band mass-to-light ratio of M/Lj 
3.6. 



^ This value is adopted 

|http://mips.as.arizona.edu/~cnaw/sun.html| 



from 
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Fig. 6. — Same as Figure[2l now including /* measurements from Gz07 (left panel) and Gi 09 ( right panel). The solid light grey region in 
the left panel shows our the systematic error margin for assuming a Chabrier IMF (see S I4.3|l . The hashed region shows the systematic 
error margin for /* assuming a Salpeter IMF. Yellow diamonds in the left hand panel represent data-points from Gz07 for the satellite 
contribution to This component can be directly compared to our low redshift satellite term (yellow dash-dash line). The total stellar 
mass fraction from Gz07 is shown by the blue squares {f^'^^+ff'^'^+fl'^^). Note that our predictions do not include the ICL contribution: 
our blue solid line corresponds to /j^'+Z^'^'"'- The right panel shows the total galaxy stellar mass fraction (/r*+/?'^°: red solid line) 
measured by Gi09. The Gi09 fit is averaged over all groups at z < 1 and so is most comparable to our Z2 bin (solid blue line). In both 
cases, our results suggest significantly lower values for than reported by either of these previous studies. 




Fig. 7. — Left: j+-band mass-to-light ratios of low-z group members as a function of galaxy color as inferred from our stellar mass 
estimates (Chabrier IMF). The mass-to-light ratio used by Gz07 is shown by the hashed magenta region (there is a spread in values due to 
to conversion from Cousins I to Subaru i+). Assuming a Chabrier IMF, the stellar masses inferred by using the Gz07 AI/L^+ are larger 
than ours by a factor of 2 to 10. The difference is particularly large for active (blue dash-dot line) and intermediate type (green dashed line) 
galaxies. Right: i+-band mass-to-light ratios assuming stellar masses are increased by 0.25 dex. To first order, these are the mass-to-light 
ratios expected for a Salpeter IMF. The difference is reduced by assuming a Salpeter IMF but the Gz07 M/L ratio is still well above the 
mean M/L ratio of our group galaxy sample. 



who find that the quenched fractions in groups is of or- 
der 40% to 60% at z < 0.5). Note however that Figure [7] 
does not include the corrections described in Section [32] 
for de-projection or the completeness and purity of the 
membership selection. These corrections were derived as 
adjustments to the total from mock catalogs that do 
not distinguish quiescent galaxies from active ones, so we 
cannot directly apply them to the separate galaxy types 



in Figure[71 To address the possibility of color-dependent 
contamination of the group sample, we can test the pu- 
rity and completeness of the membership selection using 
a subsample of galaxies with spectroscopic redshifts as 
done in Gell. The COSMOS spectroscopic redshift sam- 
ple, however, is dominated by zCOSMOS "bright" galax- 
ies with FSlAW < 22.5 and is therefore not representa- 
tive of our group membership sample {FSIAW < 24.2). 
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This will induce a degree of uncertainty that can only 
be reduced with deeper and more representative spectro- 
scopic data and/or more representative mock catalogs. 
Having noted this caveat, we use this approach to com- 
pute the completeness and purity of the membership se- 
lection within -R500C and as a function of galaxy type. 
We find that the group populations of red, green, and 
blue galaxies are overestimated by 14%, 20%, and 30% 
respectively. After applying these correction factors to 
groups in bins A and B (and for galaxies above our stel- 
lar mass completeness limit), we find that the mean red 
fraction at z < 0.48 is roughly 50%. So while we do see 
some added contamination of our group member sam- 
ple among blue galaxies, they still represent a significant 
part of the population. Studies that aim to compute fi, 
must take this varied group population into account. It is 
clear that assuming that all group members are quiescent 
could lead to large biases in estimates. 

Finally, we now investigate the impact of using a sin- 
gle mass-to-light value similar to the one use by Gz07. 
For this exercice we adopt a single value of M/Li+ = 3.8 
and we neglect the small color dependance of the M/L 
conversion from Cousins I to Subaru (see Figure [7]). 
Using the COSMOS group catalog, we re-compute /^^' 
using the absolute i+-band AB magnitude provided by 
the COSMOS photoz catalog and a mass-to-light ratio 
of M/Li+ = 3.8 to derive stellar masses. Here, we only 
consider /^^'. Because our analysis does not include ICL, 
we cannot compare with the Gz07 f^'^'^+fl'^^ compo- 
nent (these were measured together as a single compo- 
nent by Gz07). However, as noted in Figure HI most of 
the discrepancy between our results and Gz07 occurs for 
ysat^ Figure m shows that we are able to qualitatively re- 
produce the Gz07 fl'^^ using our data and M /Li+ = 3.8. 
Note however that the halo mass range and redshift range 
of our sample does not exactly match Gz07 and so Figure 
[5] is not exactly a one-to-one comparison. However, the 
fact that we obtain similar fl^^ values as Gz07 strongly 
suggests that the primary cause of the discrepancy is due 
the fact that Af/Lj+ =3.8 largely overestimates the stel- 
lar masses of group members compared to our full SED 
fitting technique. 

4.5.2. Comparison with lGiodini et oil (20031) 

Gi09 have used the COSMOS data and an earlier ver- 
sion of our X-ray catalog to derive f^,. Since we use the 
same data set and a similar group catalog, our results 
are directly comparable. The right panel of Figure [5] 
shows the comparison of our results with Gi09. Since 
Gi09 report the average over all groups at z < 1, their 
results are most comparable with our Z2 redshift bin. To 
estimate stellar masses, Gi09 have used a galaxy-type 
dependent stellar mass-to- A's band lumi nosity ratio de- 
rived by using the analytical relation from lArnouts et ahl 
assuming a Salpete r IMF. However, as p ointed 
out by lllbert et al.l (|2010[) , the lArnouts et all ()2007f ) 
M/Lk, relation was only calibrated for massive galax- 
ies. In practice the K-band mass-to-lig ht ratio varies 
with galaxy age and color. Thus the lArnouts et al.l 
(|2007f) M/Lk, relati on overestimates t he stellar masses 
of l ow mass galaxies, lllbert et al.l (|2010f) have shown that 
the lArnouts et al.l (|2007[ ) mass estimates contain galaxy- 
color dependent biases compared to SED fittin g tech- 
niques (see their Figure 28). For example, illbert et al.l 
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Fig. 8. — Impact of mass-to-light ratio assumptions on the frac- 
tion of stars in satellites, fl'^^. In addition to the default values 
from our group catalog (black error bars), values obtained from 
our group catalog using a similar i+-band mass-to-light ratio as 
Gz07 (Af/L^+ = 3.8) are shown (yellow triangles). In detail the 
two samples are not directly comparable because our sample has 
lower halo masses and is at higher redshifts than Gz07. However, 
the fact that the Gz07 results are qualitatively reproduced using 
M/L^+ = 3.8 suggests that the derivation of stellar masses is the 
primary cause of the observed difference between our results and 
Gz07. 

(|2010D find that the lArnouts et al.l (|2007l ) masses are bi- 
ased high by more than 0.3 dex for star forming galaxies 
with logjQ(M*/MQ) < 9.5. Thus, including these galax- 
ies will lead to an overestimate of 
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Fig. 9. — Impact of stellar mass estimates on In addition to 
the default values from our group catalog (black e rror bars), val- 
ues ob tained from our group catalog using the same lXrnouts et al.l 
1)20071 ) stellar mass estimates as Gi09 are shown (blue triangles). 
These are in qualitative agreement with the Gi09 results (red lines) , 
indicating that the assumptions about stellar masses are the pri- 
mary source of the difference between the two sets of results. 

Since the Gi09 results are based on COSMOS data, we 
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can make a direct comparison between their mass esti- 
mates and ours. We match our galaxy catalog with the 
photoz catalog that was used by Gi09. We then recom- 
pute /* from our catalog u sing the same proced ure as in 
Section 14.21 but using the lArnouts et all (|2007[ ) masses 
and the results are shown in Section 14.21 For this exer- 
cise we have used the same redshift and halo mass bins 
as in Section |4?2] (see Table [1]). Note however that we 
do not apply the purity /completeness and de-projection 
correction factors here. Indeed, since the Gi09 stellar 
masses are different than ours, we would need a new set 
of mock catalogs to perform this correction. The correc- 
tion values however are typically only of order 5 to 14% 
and are much smaller than difference between our results 
and Gi09 (a factor of ■~2.5). Neglecting these correction 
factors here should therefore not change our main con- 
clusions. 

Figure |9] shows that we obtain roughly similar yalues 
for as Gi09 when we use the lArnouts et al.l ()2007f ) 
masses. In detail, however, we do not exactly find the 
same mean values as a function of halo mass as Gi09 
but our methods are sufficiently different (Gell use a 
Bayesian probabilistic method whereas Gi09 use a back- 
ground subtraction method) that this is not too unex- 
pected. Part of the observed difference in Figure IH] is due 
to the fact that Gi09 use a Salpeter IMF whereas we use 
a Chabrier IMF (this accounts for about 0.25 dex). The 
remaining difference (about 0.15 dex) is probably due to 
a combination of differences in the techniques that we 
have used and to differences in stellar mass estimates. 

5. DISCUSSION AND CONCLUSIONS 

In this paper we have derived the total stellar mass 
fraction, f^,, as a function of host halo mass from z = 
0.2 to 2 = 1 using HOD methods, abundance matching 
methods, and direct estimates from group catalogs. Our 
stellar masses are derived from full SED fitting to multi- 
band photometry (including K-band) from the COSMOS 
survey. Assuming a Chabrier IMF, we find significantly 
lower estimates for at group scales than previous work. 
Including (non-IMF) systematic errors on stellar masses, 
our analysis suggests that 0.012 < /* < 0.025 at Msooc = 
10" Mq and 0.0057 < /* < 0.015 at Mgooc = 10^^ Mq. 
We will make files available upon request including our 
fi, estimates and the associated systematic uncertainties. 
Our main results are as follows: 

1. Assuming a Chabrier IMF, we find that previously 
published estimates of /* on group scales could be 
overestimated by a factor of two to five. This dis- 
crepancy is only partially reduced by assuming a 
Salpeter IMF. We investigate the cause of this dis- 
crepancy and find that a large fraction of the ob- 
served difference can be explained by the use of 
over-simplistic mass-to-light ratio estimates. 

2. We show that galaxies in groups are a mixed pop- 
ulation of quiescent, intermediate, and star form- 
ing galaxies. As a consequence, galaxies in groups 
are not well represented by any single M/L ratio 
value. The assumption that all galaxies in groups 
are quiescent, for example, will lead to biases in fi, 
estimates. 



3. We quantify the systematic uncertainty on f^, us- 
ing abundance matching methods and show that 
the statistical uncertainty on f^, 10%) is cur- 
rently dwarfed by systematic uncertainties associ- 
ated with stellar mass measurements ('~ 45% ex- 
cluding IMF uncertainties). We provide first esti- 
mates for the systematic error margin on f^, due to 
uncertainties in stellar mass measurements. While 
we have tried to provide a conservative estimate 
for these systematic errors, it is also clear that this 
is a non-trivial task that will require further inves- 
tigation. One aspect in particular that could be 
improved on in this respect would be to consider a 
larger set of SMFs that cover a more representative 
range of SPS models than we have used here. 

4. While direct measurements of using group and 
cluster catalogs are necessary and important, HOD 
and abundance matching methods can probe f^, 
over a much wider halo mass range than possi- 
ble using group and cluster catalogs. For example, 
in COSMOS, HOD methods allows us to probe 
down to the stellar mass completeness limit of the 
survey (M, - 7 x 10^ M© and Mh ~ lO" Mq at 
z = 0.37). In COSMOS, we show that the HOD 
method and the group catalog method are in good 
overall agreement but exhibit some small but in- 
teresting differences regarding the contribution to 
f^, from satellite galaxies (Z^***). 

5. Using idei itical data sets, we find that the HOD 
method of iLeauthaud et all (l2011bD and the abun- 
dance matching method of iBehroozi et all (|2010l ) 
agree to within 15% in terms of determining f^,. 
Identifying the cause of remaining small differences 
between the two methods will require further inves- 
tigation, but these small differences do not affect 
the main conclusions of this paper. 

Our finding that the total stellar fraction in groups 
is lower than previously estimated has interesting im- 
plications for the thermodynamic history of the intra- 
group gas. In simulations, the relative proportions of gas 
mass to stellar mass are directly related to the efficiency 
of cooling, star formation, and AGN feedback (e.g., 
Kravtsov et al. 2005; 'Puch wein et al.l 120081: iBode et alJ 
2009; Puchwein et al. 2010). Suppressed gas mass frac- 
tions at small radii might indicate that gas has been re- 
moved from the centers of groups and clusters by non 
gravi tational processes . For e xample, recent simulations 
from iMcCarthv et al.l ()2011[ ) have suggested that low 
entropy gas might be ejected from the progenitors of 
present-day groups by AGN feedback at high redshift 
(primarily 2 < z < 4) at the epoch when super-massive 
black holes are in quasar mode. On the other hand, low 
gas mass fractions might also be exactly compensated for 
by /*+/icL in such a way that the total baryon fraction 
is constant with halo mass and with halo radius. Dis- 
criminating between these two scenarios is key towards 
understanding how much gas may have been ejected from 
these systems. 

In Figure 1101 we have compared our results with th e 
mean gas mass fraction (/gas) from 'Pratt et al.' (20091), 
computed from a compilation of data froni .Arnaud et aLl 
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Fig. 10. — Comparison of our f^, estimates with the mean gas mass fraction from lPratt et afl 120091) . Even assuming a Salpeter IMF, our 
measurements imply that the total stellar content of dark matter halos is always lower than the gas mass fraction for halos above IQ-'^'^A^'q 
and for r < R^oOc- 
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Fig. 11. — Sum of the stellar fraction and gas fraction /gas, where /gas is the mean gas mass fraction taken from lPratt et aO 120091 ). 
The solid grey region shows the systematic error margin for /*+/gas assuming a Chabrier IMF. The grey hashed region shows the systematic 
error margin for /*+/gas assuming a Salpeter IMF. 



(pool . IVikhlinin et all (l2006l) . and ISun et al.l (l2009h . 
This result demonstrates that /* is always lower than 
/gas for halos above 10^^ Mq, for the assumption of either 
a Salpeter or Chabrier IMF. Figure [TT] shows the sum of 
and /gas • Here the solid grey region indicates the sys- 
tematic error margin for /^+/gas assuming a Chabrier 
IMF; the grey hashed region shows the systematic error 
margin for /*-l-/gas assuming a Salpeter IMF. Even in 
the latter case, we find a shortfall in /*+/gas compared 



to the cosmic mean, which increases towards lower halos 
masses. 

It is important to note that Figure [TT] does not include 
a baryonic contribution from ICL. While including ICL 
in the baryon budget is obviously an important issue to 
be adressed, /icl estimates suffer from large uncertain- 
ties related to the choic e of a M/L value. Fo r exam- 
ple, the ICL estimates of iGonzalez et all (|2007l ) assume 
M/Lj = 3.6; Figure [7] demonstrates this is on the high 
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end compared to MjLi values predicted from SPS mod- 
els. Finally, we note that the main claim in this paper, 
that ]i, estimates are lower than previously estimated, is 
dominated by the /^^' component of Z^. 

In summary, we have presented new measurements of 
the stellar fraction /j, in groups and clusters from sev- 
eral complementary approaches, which we believe to be 
more robust than previous estimates. However, precise 
measurements of are limited by our systematic un- 
certainties associated with stellar mass measurements. 
Future improvements to the /j, measurements presented 
in this paper will benefit most from efforts that lead to 
an improved understanding of galaxy stellar masses. 
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